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As the Einstein equations are non-linear, spatial averaging and temporal evolution do not com- 
mute. Therefore, the evolution of the averaged universe is affected by inhomogeneities. It is, 
however, highly controversial how large these cosmological backreaction effects are. We use the 
supernova data of the Constitution set up to a redshift of 0. 1 in order to analyse to what extent 
the measurement of the Hubble constant is affected. The size of the effect depends on the size of 
the volume that is averaged over. The observational results are then compared to the theory of the 
backreaction mechanism. 
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1. Introduction 

In cosmology it is very common to assume a homogeneous and isotropic universe. As we 
know that there exist structures in the universe [p, this homogeneity and isotropy can only be 
statistical on large scales but not exact. Due to the non-linearity of Einstein's equations spatial 
averaging and temporal evolution do not commute. Thus, local inhomogeneities can affect the 
expansion of the background universe via the so-called backreaction mechanism |g[ ffl, g, gj. In 
this work, we probe the influence of backreaction on the measurement of the Hubble rate using 
supernova type la data. 



2. Averaging 

Many observables are averaged quantities. As all observed objects lie on our past light cone, 
it would be appropriate to average the observables over this light cone. This is however a very 
difficult task that has not been achieved yet. Instead one can use spatial averages at the cost of 
introducing an error. In order to keep this error at an acceptable level, it is necessary to limit the 
use of spatial averages to low redshifts. In our analysis we use supernovae up to z = 0.1. In that 
range, this averaging method is justified. 

The averages of observables are calculated within a certain domain D. Its volume is given by 

V D (t) = J D W D {x)y/d^dx, (2.1) 

where Wd(x) is the window function specifying the domain. Then the spatial average of an observ- 
able O within D is 

(0) D = ^f D W D (x)0(t,x)^gTjdx. (2.2) 
An effective scale factor ap can be defined via the domain volume: 

where the subscript denotes the present time. The effective Hubble rate then determined by 

H D = ±. (2.4) 
a D 

Following Buchert's formalism, the Einstein equations can be averaged in order to obtain the 
effective Friedmann equations for a dust universe [j|] : 

d D \ 2 87TG 

Peff , (2-5) 

(Peff + 3peff) • (2.6) 



O-DJ 3 



an 3 

These equation include the energy density and pressure of an effective fluid, which are given by 

Pas = {p)d--^({Q)d+{S)d) , (2-7) 
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(Q)d denotes the kinematical backreaction and {M)d the averaged spatial curvature. 

We want to analyse the influence of backreaction effects on the measurement of the Hubble 
rate. Its average value obtained by observing objects within a domain D is denoted as Hp. Assum- 
ing that there exists a global value Ho, we can define the fluctuation of the Hubble rate as 

_ H D - H 

oh = — • (2.9) 

^0 

Without backreaction the average value <5# equals zero. Considering backreaction effects, the value 
of 8h becomes slightly negative for small domain sizes. But the main effect is that the variance of 
Sh is increased compared to the case when backreaction effects are not taken into account. 

3. Gaussian window function 
3.1 Method 

It is not obvious which choice of window function yields the best results for a test of backre- 
action effects. The first try was to assume a spherically symmetric domain described by a gaussian 
window function 

W D (r) = ^— expf-— T ) , (3.1) 



where Rp specifies the size of the domain. Then the variance of Sh due to backreaction effects is 
given by [@]: 

where Rh is the Hubble radius. The values for the power spectrum 9*^ are taken from WMAP5 
measurements [Q]. 

For the analysis, we use supernova type la data from the Constitution set [P] up to a redshift 



of 0.1. We used the data that were fitted with SALT2 [10]. The number density of SNe needs to 
be approximately constant within the considered domain. Thus, the number of SNe N(r) in the 
distance interval [r, r + dr] must be proportional to r 2 Wo(r). That means that we have to choose a 
subset of SNe, whose distribution in space corresponds to that of the considered window function. 

Figure [jj shows the distribution of the 178 SNe in the Constitution set up to redshift 0.1. 
For the test of backreaction effects, it is essential to compare the value of 8h of differently sized 



domains. The domain size can be changed by varying Ro as given in equation ( J3- ![ )- We chose 
to use five domains. Their corresponding window functions r 2 Wo(r) are plotted in the figure, the 
black curve being the sum of these functions. However, the window functions only determine how 
many SNe at a certain distance are assigned to a subset, but not which individual SNe. The actual 
assignment of SNe to a subset within one realisation is done randomly, but in a way that all subsets 
are disjoint and thus statistically independent. As soon as one uses many realisations, i.e. different 
assignments of SNe to the subsets, the statistical independence is lost. 

In order to calculate the Hubble rate Hd, we need to know the redshift z and the distance 
modulus /i for each supernova. The distance modulus, however, depends on the calibration of the 
absolute magnitude of the SNe. Thus, a different calibration leads to different values of 8h, if one 
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Figure 1 : Gaussian window function, (a) Distribution of the SNe and r 2 Wp (r) for five different val- 
ues of Ro = 45,60, 80, 100, 120. The black curve is the sum of all five window functions r 2 Wo(?"). 
(b) 8h obtained from ten different realisations of subsets. Also shown are the variances with and 
without backreaction. 



uses a global Hq that was not obtained by using the same data set, but by other observations such 
as WMAP. It turned out that the test result is very sensitive to the calibration. Therefore, we have 
to determine Hq using the same data set as for calculating Hp. 



3.2 Results 



Figure lb shows 5# for ten random realisations of subsets. The five data points obtained from 
the five subsets in each realisation are connected by lines, r is the average distance of the SNe in 
a subset. The pink curves indicate purely the measurement errors of the SNe. For the red curves 



-y/Var (5#) from equation (3.2) is added in quadrature to the measurement errors. Here, the domain 



scale Rd has to be expressed in terms of the average distance as Rd = ^J%j%r. The global Hubble 
rate that is needed for the calculation of 8h was chosen such that the data points at large distances 
lie within the variance limits. There is a trend of increasing 8h with decreasing distance. The model 
with backreaction effects seems to describe the data better than the one without these effects. 

In order to test backreaction effect, we need however a more quantitative analysis. Therefore, 
we determined the optimal Hq for each realisation once with and once without backreaction effects. 
Then we calculated the likelihoods of each model given the data. In 27 out of 100 realisations the 
model including backreaction effects was favoured. However, in none of the realisations one of the 
two models was favoured significantly. 



4. Tophat window function 

As it was not possible to detect backreaction effects using a gaussian window function, we 
tried a different ansatz. Since backreaction effects are larger at smaller distances, our aim was to 
minimize the distance of the first data point. This can be achieved by only using the nearest SNe to 
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Figure 2: Tophat window function, (a) Distribution of SNe. (b) 8h obtained from the SN subsets. 
Also shown are the variances with and without backreaction. 



calculate the first data point. So we binned the SN data according to distance, where the binwidth 
is increased with increasing distance (see figure ^a|). The corresponding window function is 

W D (r) = e{R D )@ (t-^Rd) • (4-1) 
Then the variance of <5# is given by 

Val( « = 153664^ (iTF? {£) L **/*>VfcM* ■ <«> 

The relation between Rp and the average distance is Rp = 49/68 r. 



An advantage of the new window function ( |4.1| ) is that the assignment of SNe to subsets 
is unique. So we do not need to consider different realisations. Like in the previous case, we 
determined the optimal global Hubble rate Hq for the model with and that without backreaction 
and subsequently calculated 8h for each subset. The result is shown in figure Note that the data 
points for the two models differ slightly as we have used different values of Hq. The quantitative 
analysis shows that the model with backreaction is only about twice as likely as the models without 
backreaction effects. Thus, we have not found any evidence for backreaction. 



5. Conclusion 



Theoretically, backreaction influences the measurements of the Hubble rate by increasing its 
variance. This effect should be observed, if it was possible to measure the Hubble rate at different 
locations in the universe. We are, however, restricted to our local universe. So it is possible that 
our local measurements are by chance consistent with a model that does not include backreaction 
effects. If that was the case, we would not be able to detect those effects using the test presented 
in this work. Therefore, the test can potentially prove the existence of backreaction effects, but it 
cannot prove that there are no such effects. 
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Using the currently available supernova data, we could find some slight hint of backreaction, 
but no evidence. Larger data could help providing that evidence. A larger number of supernovae 
leads to a smaller variance in the model without backreaction. If the measured values of 8h stayed 
approximately the same for a sufficiently large data set, then the data would lie significantly outside 
the variance limits of a model without backreaction. In that way future data sets have the potential 
of providing the evidence for the backreaction mechanism. 
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